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Abstract 

In the following some aspects of inclusive hard processes in photon induced reactions are 
reviewed. After a discussion on the properties of hard processes, the phenomenology of jet 
production and of charmonium production is presented in the context of the next-to-leading 
logarithm approximation of QCD. 
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1 Introduction 


Over the last few years the physics of hard processes in photon-induced reactions {e.g. photo¬ 
production at HERA and photon-photon collisions at TRISTAN and LEP) has almost reached 
the status of sophistication of hard processes in purely hadronic reactions. The experimental 
results are becoming more and more accurate and many observables measured in hadronic 
reactions are now being accessible in photon-induced reactions while the corresponding theo¬ 
retical calculations have been performed at the next-to-leading order (NLO) of perturbative 
QCD. Since PHOTON ’95 [1J an enormous progress has been made in this respect so that 
one should be able to start quantitative phenomenology. The main points that deserve special 
attention now are: 1) matching the calculated observables with the experimental ones; 2) 
matching the appropriate non-perturbative input to the conventions used in the NLO calcu¬ 
lations. 

It is well known that the physics of photon-induced processes is more complex than that 
of pure hadronic reactions. This arises from the fact that the photon acts either as a parton 
which couples directly to the hard scattering and contributes to direct processes, or as a 
composite object (a bag of partons) the constituents of which couple to the hard scattering 
via the photon structure functions: this leads to resolved processes. Furthermore, as discussed 
in ref. 0 the photon structure function contains two pieces: 

F%h{x,Q) = Ffff(x,Q) + i = quark, gluon (1) 

where the first term on the right hand-side is similar to the hadronic structure function 
while the second term increases logarithmically with Q 2 and is asymptotically calculable in 
perturbation theory. 

In the following, in a theoretical introduction, the structure of a “typical” cross section 
will be derived stressing the features distinguishing the direct processes from the resolved 
ones. Two topics are selected for further discussion: the production of jets and single particle 
on the one hand and the production of hidden charm for which HERA offers the possibility 
to test models of J/T production at the Tevatron. 

2 Photon induced reactions: theory 

We consider photon-photon collisions as an example. The lowest order diagram is shown in 
fig- 0a- To apply the perturbative approach a large scale is needed which is provided by the 
transverse momentum of the produced jet or the mass of the heavy quark. Considering jet 
production, the Born approximation is valid for large values of p T /\fs where unfortunately the 
cross section is very small since it behaves as do^ et jdpT ~ a 2 /p\.. To obtain reliable predictions 
at lower p T values it is necessary to consider diagrams with an extra gluon emission, leading 
in principle, to 0(a s ) corrections to the cross section (fig. |b, c). When integrating over the 
phase-space of final state partons to reconstruct e.g. the single inclusive jet cross section, one 
encounters “dangerous” regions where the virtualities of some fermion propagators may vanish 
and the corresponding matrix elements are not defined. This happens, for example, when the 
final state parton of momentum p l is collinear to the initial photon of momentum k a (s ai = 
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Figure 1: Feynman diagrams for 77 scattering, (a): Born approximation; (b) real corrections; 
(c) virtual corrections. 


(k a — pf) 2 —»► 0 ) or the emitted gluon is collinear to a final state quark (s 13 = (p, + p 3 ) 2 —>► 0 ). 
Considering the s ai singular region, the process can be pictured as in fig. ^]a where the initial 
photon fragments into a collinear qq pair followed by the hard 2 body —> 2 body scattering 
of the q or q, the other parton flying off down the beam pipe. The photon fragmentation is a 
soft process (the relevant scale s ai is small) which takes place on a long time scale before the 
short-distance process scatters the partons at large transverse momenta. The cross section 
can be written as 


da™ 

dpxdg 


dzP q / 1 (z) In 



da™ 

dprdr] 


+ k d ( Pt ,m ), 


( 2 ) 


where A is a cut-off introduced to regularise the collinear singularity, and M ~ 0(p T ) is an 
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Figure 2: Diagrammatic representation of the collinear singularities. 


arbitrary scale introduced to separate the contribution from the collinear region, generating 
the “large” ln(M 2 /A 2 ) piece, from the hard region (rouhgly s ai > M 2 ) leading to the term 
K D containing only finite terms such as ln(p T /M 2 ). The above expression is rigourously 
independent of the scale M. I 11 fact, at this level of approximation the scale M is not necessary 
but it is introduced for later use. The integration variable in eq. is simply the fraction of 
the photon momentum carried by the interacting q or q. Using the factorisation theorem, we 
can substitute to the “large logarithm” term P q /^(z) ln(M 2 /A 2 ) the photon structure function 
F q /j(z,M), measured in 7*7 collisions, which satisfies evolution equations of type: 

= P L1 + y p tj ® Fjh (M), i,j = q,g (3) 

where the symbol <g) denotes a convolution over the longitudinal variable and the Pij are 
the usual splitting functions. After replacing P g / 7 ln(M 2 /A 2 ) by F q / 1 (M) in eq. (||) the 
compensation in the M dependence on the right-hand side is only approximate as the evolution 
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equation effectively resums large higher order corrections (of type [a s (M) ln(M/A QCD )] n ) in 
F q /y(M) whereas no such resummation is performed in K D . The cross section can be written 


da 77 
dpTdrj 


( da 11 (0) 

\dpxdi] 


Os 

2n 


K d ( Pt ,M )J + J dzF qh (z,M) 


da 79 
d'prdfi) 


( 4 ) 


where the K D term is a correction to the Born cross section da 11 ^ and contributes to the 
direct cross section evaluated in the NLO approximation. The last term in the equation is the 
resolved process. Since in the latter case only part of the energy of the photon contributes to 
the large p T process, the rest being carried by the longitudinal fragment, the resolved term 
is qualitatively (be. in the LO approximation) distinguished from the direct one by: 

- a softer p T spectrum; 

- a jet system boosted in the direction opposite to that of the resolved photon: e.g. a backward 
moving resolved photon produces forward going jets; 

- some hadronic activity in the direction of the resolved photon. 

In the NLO approximation the direct and resolved components are not separatly observable 
since they depend on the arbitrary renormalisation scale M. However, using the properties 
above it is possible to define observables related to these two components (see below). 

Until now we have discussed features related to the initial state singularities. The final state 
singularities (of type s 13 —> 0 ) also lead to 2 body —> 2 body hard scattering with the singular 
behaviour associated to the fragmentation process exactly as in purely hadronic processes (fig. 
|b). The cross section is made finite by adding the virtual diagrams and properly defining 
the jets, merging the two almost collincar partons into one jet, or by convoluting with a 
fragmentation function and using the factorisation theorem to build the scaling violations in 
the fragmentation function. 

One could continue the perturbative analysis and consider the emission of more partons: 
this will lead to double resolved processes where each photon interacts via its structure function. 
Finally, the 77 cross section for inclusive jet production takes the form: 

da 11 _ da D | da SR | da DR 
dflrdr] dpTdrj dfirdr) dpTdr} 

with the direct , single resolved and double resolved pieces given in the NLO approximation 
by an expansion of type: 


da D 

dpxdp 


(R) 


da 11 (0) 
dpxdg 


O^s 

27T 


K u (R;p t , M), 


( 6 ) 


and similarly for the other terms. The variable R specifies the jet cone size. The factori¬ 
sation scale variation associated to the inhomogeneous term P n in eq. (|]) is compensated 
between the D, SR and DR pieces while the variation associated to the homogeneous terms 
Pij compensates within each of these pieces between the lowest order and the correction terms. 
Thus only eq. @ is expected to be stable under variation of scale M but not its individual 
components. 

For hadron production an extra convolution of the above expresssions with the relevant frag¬ 
mentation function is needed. For photo-production processes one of the F l /„,(x, M ) should be 
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replaced by the hadronic structure function. In an obvious notation, the cross section takes 
the simpler form: 


da lp da D da R 

dpTdij dprd/rj dpTdr) 


( 7 ) 


3 Phenomenology of jet production 


One of the aims of studies of hard processes in 7 induced reactions is the NLO determination of 
the photon structure function ||. To achieve this it is very useful to isolate in an experimental 
way the resolved component of the cross section. In a direct process, all the photon energy 
is given to jets at large p T so that by momentum conservation one has for a 2 -jet event (the 
photon is assumed to move towards negative rapidity): 


= 


E^e T E^e 
2EL 


= 1 


( 8 ) 


Higher order correction are not expected to change this relation drastically. In contrast, for a 
resolved event x 7 < 1 since part of the photon energy disappears in the beam pipe (fig. @a) (in 
ref. Q a modified definition of x 7 is proposed). At HERA, the experimental groups choose the 
value x 7 = .75 as a cut-off so that above this value the events are produced mostly by the direct 
process while below they are mainly resolved. The integrated di-jet cross sections for HERA, 




Figure 3: Integrated di-jet cross section in the domains 0 < x 7 < 1 (solid line), x 1 < .75 
(dashed line), x 7 > .75 (dotted line). From ref. & 

shown in fig. nicely illustrate, in the NLO approximation || the different characteristics 
of the two components. In the first figure the direct component shows a faster decrease with 
Et, while in the second one the resolved cross section is seen to contribute mainly forward 
jets. In fig. [I] a very good agreement is seen, between the NLO theory |>j and experiment |6j], 
in the shape of the di-jet angular distribution (the angle is measured in the di-jet rest frame): 
as expected the resolved component has a steeper angular dependence due to the importance 
of gluon exchange diagrams while for the direct term only quark exchange is allowed. In the 
next figure, we compare theoretical predictions of two independent groups 1,0 with ZEUS 
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Figure 4: Dijet angular distribution as measured by ZEUS normalized to one at cos 6* = 0 
compared with LO (dash lines) and NLO result (solid lines), from ref. SI- 

data for the di-jet cross section || da/dfj, integrated over the phase space E Tl , E T2 > E 0 and 
\ ,} li — 1 < -5, with i] = ( r/i + 77 2 )/2. The first point to notice is the disagreement between the 

two sets of theoretical predictions: this arises because the groups do not calculate the same 
observable due to instabilities in the perturbatively calculated quantity. This is related to the 
fact that the boundary condition Et x = Et 2 = Eq is an infrared singular point. Typically, in 
the NLO approximation one needs to consider the production of 3 partons in the final state: 
it can be separated into two classes: 

a 2 - 3 = a^(y c ) + a^\y c ) (9) 

with the parameter y c such that the first term on the right hand side contains events generated 
with “dangerous” invariant masses [i. e. leading to soft and collinear singularities in the matrix 
element] Sij > y c s ab , while crjj^ 2 (y c ) contains events which look like 2-body hard scattering as 
discussed above. This part contains all divergences and y c is chosen small enough ( y c = ICO 5 
to ICO 2 ) so that approximations can be used to extract divergences analytically. Upon adding 
the virtual corrections which are also of type 2 body —> 2 body scattering all divergences 
cancel and the theoretical cross section is: 

a = a^(y c ) + a 2 ^(y c ) ( 10 ) 

each piece being regular, the first one beeing calculated numerically and the second one semi- 
analytically. By histogramming one reconstructs any observable but it should be checked that 
the result is independent of y c . ft turns out that for the dijet observable of fig. [|, the condi¬ 
tion Et x , Et 2 > Eq introduces constraints on the phase space which spoil the y c compensation 
between the terms in the above equation. Harris and Owens || (who use a more elaborate 
method than the one described above) observe that the remaining dependence on y c is much 
less than the experimental error bars and present results with some value of y c while Klasen 
and Kramer 0 modify the boundary condition to E Tl > E 0 , allowing for a smaller E? 2 if the 
third unobserved jet has transverse energy less than Et 3 < 1 GeV. This is sufficient to remove 
the y c dependence but the result is rather sensitive to the latter energy cut. As seen in the 
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Figure 5: Dijet cross section vs. fj integrated over E ] fi > E 0 for E 0 = 6, 8,11,15 GeV; the 
four figs, on the left are as for x 1 > 0.75; those on the right are for x 1 < 0.75. The data are 
from /f|/. Thin lines from ref. m thick dashed lines from ref. 07 - 


figure this slight modification of the boundary condition considerably affects the predictions 
||. If one considers that the calculation of ref. [!5] is closer to the experimental observable 
there remains a drastic disagreement betwen the data and the theoretical predictions specially 
concerning the resolved cross section. Disagreement between theory ll-i and experiment 
|T^j is also seen when comparing the rapidity distribution of single inclusive jet production 
as displayed in fig. |6|: this is particularly marked at low E T and large pseudo-rapidities, i.e. 
where the resolved component plays a dominant role. This leads us to the important question 
of whether or not the theoretical predictions are able to match the experimental observables. 
Two questions can be raised: 

- is it possible to match theoretical jets, made up of partons, with experimental jets recon¬ 
structed, via various algorithms, from energy deposited in calorimeter cells? 

- how does one take into account the transverse energy of the underlying event, i. e. the energy 
generated by the interaction of the spectator partons? 

The first point was studied many years ago in connection to jet production at the Tevatron 
and was recently re-analysed in detail for HERA by Kramer and collaborators Jl4|: the 
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idea is to introduce an extra parameter R sep controling the width of the jets which is adjusted 
to fit the data: a better description of the data by the NLO calculations is indeed achieved at 
the expense, however, of the predictive power of the theory. 

Concerning the second point the problem can be qualitatively understood in the following 
way. At HERA, in events with direct photo-production of jets the transverse energy of the 
underlying event is produced by string-like effects: typically a transverse energy of 300 — 400 
MeV per unit of rapidity is expected. On the contrary, in resolved photo-production, the rem- 
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Figure 6 : Comparison of the single jet inclusive cross section measured by HI JZ|/ compared 
with NLO predictions ||/ using 2 different sets of photon structure functions. Transverse 
energy of the underlying event at HERA as measured by HI \ H\] , 


nants of the photon can undergo soft or semi-hard interactions with those of the proton |L5 


generating transverse energies in the GeV range: the larger the energy of the photon remnant 
(he. the smaller is) the larger is the underlying transverse energy: this is illustrated in 
fig- | from HI ]n| which shows the Et outside the jets as a function of x 7 . Such an effect 
is of course present in purely hadronic reactions but is relatively less important because the 
much higher Et values of jets probed in hadronic colliders. On the other hand, jets in 77 
reactions should be less affected except at very low E T values ( E T < 5 GeV at LEP1) where 
the double-resolved process is dominant. Good agreement of TRISTAN |16!. [L7| and LEP [|l^ 
data with the NLO predictions |19| has been obtained for inclusive two-jet production. 

A discussion about the determination of the gluon density in the photon using jet cross 
sections is given, in the LO approximation, in ref. [Q. 

A way to avoid problems related to jet definition or to underlying transverse energy is 
to consider single hadron production. Good agreement is indeed found in photo-production 
between the NLO predictions |20| and HERA data pl[. 


4 Charmonium production 

There has been new developments concerning the production of hidden heavy flavor at the 
Tevatron where the usual model predictions, based on the color singlet model, for prompt T 
production fall an order of magnitude below the experimental results. The interest of HERA 
lies in the fact, that many of the new parameters introduced in the model can be tested 
independently. We consider here only the non-diffractive mechanism, corresponding to a T 
inelasticity factor of z = p.k^/p.k^ < .9. In the factorisation approach of ref. the cross 
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Figure 7: Some Feynman diagrams for charmonium production, (a): colour singlet production; 
(b): diagram contributing only to colour octet production. 


section for the production of a heavy quark bound state H is 

daQ/p —> HX ) = '^daQ/p —» QQ[n]A") < O h [ n] >, 


( 11 ) 


where da is the cross section for producing a heavy QQ sytern in state [n] defined by its color 
[1, 8], its spin [0,1] and orbital angular momentum. This term is calculable perturbatively, the 
large scale beeing provided by the heavy quark mass or eventually by the transverse momentum 
of the system. The factor < O h { n] > describes the non-perturbative transition from the 
state with quantum numbers [n] to quarkonium H. Typically < C^pSi, 1] >= 1.16 GeV 3 
(from T leptonic decay) while the 8 matrix elements < [ 3 5i, 8] >~ < 0' 3/ [ 3 S'o,8] >~< 

8 ] > /ml ~ 1CT 2 GeV 3 as determined from hts to the Tevatron data |23fl . The colour 
8 values are consistent with non-rclativistic QCD which predicts their suppression by powers 
of the velocity of the heavy quark in the bound state |24|]. Although the 8 matrix elements 
are small it may happen that the perturbative cross section dcr([8]) is large so that a non 
negligible contribution occurs. In the colour 1 model the 8 matrix elements are assumed 
to vanish. Many processes contribute to the cross section dd(QQ[n\): it is convenient to 
distinguish fusion processes which are important at small p T , but are suppressed by a factor 
m 2 /p 2 , from fragmentation processes where the QQ state is found in the decay of a Q or 
gluon jet produced at large transverse momentum. Furthermore, both processes come in the 
direct variety if the photon couples to the hard sub-scattering, or the resolved one if the photon 
structure function is involved. These mechanisms have different z and p T dependences so that 
varying the kinematical conditions they could be separated. 

Consider first direct fusion processes which should be a good approximation of cross sections 
at small transverse momentum ( p T < m Q ) or integrated over p T . The basic colour singlet 
diagram is shown in fig. [Tja : the heavy QQ pair is produced by photon-gluon fusion and an 
extra gluon is necessarily emitted so that a QQ state in a 1 component can be projected 
out. NLO corrections to this process have also been calculated [25]. A similar diagram 
contributes also to the production of a colour 8 state together with the diagram of fig. 0b: 
because of the gluon exchange this diagram will be enhanced near the kinematical region 
z ~ 1. The comparison of the theory ]26] with HERA j [27| . [28| | data is shown in fig. [| 
where one sees the excellent agreement between data and the colour 1 NLO predictions while 
the 8 component seems to yield much too large a contribution at large z values. To claim 
quantitative disagreement is premature as the non-perturbative 8 matrix elements may have 
been overestimated in the Tevatron analysis; furthermore, doubts have been raised about the 
validity of the factorisation approach (velocity expansion) in the large z region [|30] . 


The resolved fusion diagrams are easily obtained by substituting to the photon a gluon, 






















Figure 8 : Fusion contribution to charmonium production at HERA. 


fragment of the photon: obviously this process dominates at small z where it is further found 
that the 8 contribution overwhelms the 1 one by more than one order of magnitude: thus the 
region at small p T , small z or equivalently large rapidity since z ~ e~ y iab, should be able to 
probe the 8 matrix elements. 

At large p T , on the other hand, the production of charmonium in the fragments of a jet 
should be considered mM- In the factorisation approach, the fragmentation funtion of 
parton i into H takes the form: 

= 5>^ cc - [n] (z,M) < O h [ n] >, i = c,g (12) 


where the functions d, are perturbatively calculated (all scales involved are large) and the 
same non-perturbative matrix elements as above appear. A detailed study at the NLO order 
[5lJ, including both direct and resolved processes, has recently been performed where it is 
found an overwhelming dominance of the colour 8 channels at HERA for large y Uh and p T , 
the more so the larger the yp invariant mass. 

The exclusive channel yp —> Ty has also been proposed Jl3| as a test of the model as it is 
obviously dominated at large z by the 8 component (see fig. j7a with the final state gluon 
changed into a photon). 

Thus HERA appears a promising place to support or invalidate the colour 8 model if detailed 
measurements of specific channels can be done over a wide rapidity range for JT production. 
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